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The ground-state energy of the doubly magic nuclei ^He and ^'^O has been cal- 
culated within the framework of the Goldstone expansion starting from mod- 
ern nucleon-nucleon potentials. A low-momentum potential Vio-ij-h has been 
derived from the bare potential by integrating out its high-momentum compo- 
nents beyond a cutoff A. We have employed a simple criterion to relate this 
cutoff momentum to a boundary condition for the two-nucleon model space 
spanned by a harmonic-oscillator basis. Convergence of the results has been 
obtained with a limited number of oscillator quanta. 

1. Introduction 

As is well known, the strong repulsive components in the high-momentum 
regime of a realistic nucleon-nucleon (NN) potential Vnn need to be renor- 
malized in order to perform perturbative nuclear structure calculations. 
In Refs. [1,2] a new method to renormalize the NN interaction has been 
proposed, which consists in deriving an effective low-momentum potential 
Mow-k that satisfies a decoupling condition between the low- and high- 
momentum spaces. This Viow-k preserves exactly the on-shell properties of 
the original Vnn up to a cutoff momentum A, and is a smooth potential 
which can be used directly in nuclear structure calculations. 

In the past few years, we have employed this approach to calculate 
the ground-state (g.s.) properties of doubly closed-shell nuclei within the 
framework of the Goldstone expansion [3,4], using a fixed value of the cutoff 
momentum. 

Recently, we have investigated how the cutoff momentum A is related to 
the dimension of the configuration space in the coordinate representation 
[5], where our calculations are performed. We have shown how the choice 
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of a cutoff momentum corresponds to fix a boundary for tlie two-nucleon 
model space. 

In tlie present work, we calculate the g.s. energy of ^He and ^^O in the 

framework of the Goldstone expansion with different TV TV potentials. To 
verify the validity of our approach, we compare the ''He results with those 
obtained using the Faddeev-Yakubovsky (FY) method. 

The paper is organized as follows. In Sec. 2 we give a brief description 
of our calculations. Sec. 3 is devoted to the presentation and discussion of 
our results for ''He and ^^O. A summary of our study is given in Sec. 4. 



2. Outline of calculations 

As mentioned in the Introduction, the short-range repulsion of the NN 
potential is renormalized integrating out its high-momentum components 
through the so-called Mow-k approach (see Refs. [1,2]). The Vjow-k preserves 
the physics of the two-nucleon system up to the cutoff momentum A. While 
this low-momentum potential is defined in the momentum space, we per- 
form our calculations for finite nuclei in the coordinate space employing a 
truncated HO basis. This makes it desirable to map the cutoff momentum 
A, which decouples the momentum space into a low- and high-momentum 
regime, onto a boundary for the HO space [5]. 

If we consider the two-nucleon relative moticm in a HO well in the mo- 
mentum representation, then, for a given maximum relative momentum A, 
the corresponding maximum value of the energy is 



-E'max = ) (1) 

where M is the nucleon mass. 

We rewrite this relation in the relative coordinate system in terms of 
the maximum number TVmax of HO quanta: 



for a given HO parameter fuv. The above equation provides a simple cri- 
terion to map out the two-nucleon HO model space. Let us write the two- 
nucleon states as the product of HO wave functions 



\a b) = InJaja, ribhjb) 



(3) 
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Wc define OTir HO model space as spanned by those two-nucleon states that 
satisfy the constraint 



2na + la + ^rib + h < N^e^ . (4) 

In this paper, making use of the above approach, we have calculated 
the g.s. energies of ^He and ^^O within the framework of the Goldstone 
expansion [6]. We start from the intrinsic hamiltonian 



^ ^ i=l i<j 

where Vij stands for the renormalized Vnn potential plus the Coulomb 
force, and construct the Hartree-Fock (HF) basis expanding the HF single 
particle (SP) states in terms of HO wave functions. The following step is to 
sum up the Goldstone expansion including contributions up to fourth-order 
in the two-body interaction. Using Fade approximants [7,8] one may obtain 
a value to which the perturbation series should converge. In this work, we 
report results obtained using the Fade approximant [2|2], whose explicit 
expression is 



r,,|„i ^0(l + 7l+72) + ^l(l+72) + ^2 

1+71+72 



where 



^ E2Ei - Ej ^ _Ei±Ern 

'Yl El p ip2 ' '''2 T71 ' 

II1III13 — hi^ £'2 

Ei being the ith order energy contribution in the Goldstone expansion. 

Our calculations are made in a truncated model space, whose size is 
related to the values of the cutoff momentum A and the two parameter. The 
calculations are performed increasing the A^'max value (and consequently A) 
and varying fwo until the dependence on A^max (A) is minimized. 



3. Results 

We have calculated the binding energy of ^He using different Vat at 's, and 
compared our results with those obtained by means of the FY method. This 
comparison is made in order to test the reliability of our approach. In Figs. 
1, 2, and 3 the calculated **He g.s. energies obtained from the CD-Bonn [9], 
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N^LO [10], and Bonn A [11] iVTV potentials are reported, for different values 
of hhj, as a function of the maximum number -/Vmax of HO quanta. The FY 
result [12,13] is also shown. 
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Fig. 1. Ground state energy of ''He calculated with the CD-Bonn potential as function 
of Afmax, for different values of hu). The straight line represents the Faddeev-Yakubovsky 
result, while the dashed one our converged result. The difference in energy between the 
latter and the Faddeev-Yakubovsky result is also reported. 



-24.5 

_-25.0 
> 

CD 

— -25.5 

LU 

-26.0 



h\0 'He 


0-39 MeV 






* Kco=30 






▼ fea)=29.5 






• S 00=29 






■ Sco=28.5 






♦ /;(0=28 


2.0 4.0 


6.0 


8.0 




N 






max 




Fig. 2. Same as Fig 


1, but for the N^LO potential. 



For the sake of clarity, in Table 1 we report the numerical values obtained 
with the CD-Bonn potential. From the inspection of Table 1 it can be seen 
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Fig. 3. Same as Fig. 1, but for the Bonn A potential. 



that the g.s. energy does not change mcrcasing iVmax from 4 to 6 for hw ~ 36 
MeV. On these grounds, we choose as our final result that corresponding 
to the above hio value, i.e. -25.92 MeV. Moreover, we find it worthwhile to 
introduce a theoretical error due to the uncertainty when choosing hujconv, 
which corresponds to the one with the faster convergence with A^max- We 
estimate this error as the largest difference in energy between the final result 
and those corresponding to the two fku values adjacent to fkUconv For the 
CD-Bonn potential, we see that this difference is 0.05 MeV for the largest 



Table 1. Ground state energy of *He (in MeV) 
calculated with the CD-Bonn potential for differ- 
ent values of hw and A^max 



riLO Afmax = 2 Amax = 4 Amax = 6 

(in MeV) 



35.0 -26.11 

35.5 -25.10 

36.0 -25.86 

36.5 -25.71 

37.0 -25.56 



-26.13 -26.04 

-26.03 -25.97 

-25.92 -25.92 

-25.83 -25.87 

-25.71 -25.80 



Similarly, the results for the "^He g.s. energy with the N'^LO and the 
Bonn A potentials are (-25.02 ± 0.05) and (-27.78 ± 0.03) MeV, respec- 
tively. These values are in good agreement with the FY results, the largest 
discrepancy being 0.39 MeV for N'^LO potential. 
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Fig. 4. Ground state energy of ^^O calculated with the CD-Bonn potential as function of 
A^max, for different values of hu). The straight line represents the experimental value [14], 
while the dashed one our converged result. The difference in energy between the latter 
and the experimental value is also reported. 



We have also calculated the g.s. energy of ^^O starting from both the 
CD-Bonn and the Bonn A potential, as reported in Figs. 4 and 5, respec- 
tively. With the CD-Bonn potential, the converged value, obtained for hw = 
27.25 MeV, is equal to (—117 ± 1) MeV, the discrepancy with the experi- 
mental value [14] being 11 MeV. This value is slightly different (« 1 MeV) 
from the one reported in our previous paper [5], because in the present 
work we have decreased by a factor 2 the spacings between the hw values. 
It is worth to point out that this result is consistent with those obtained by 
Fujii et al. using the unitary model-operator approach [15], and by Vary et 
al. in the no-core shell model framework [16]. 

A better agreement with experiment is obtained using the weaker tensor 
force NN potential Bonn A, our "'^^O g.s. energy being (—130.0 ±0.5) MeV. 
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4. Summary 

In this work, wc have calculated the g.s. energy of the doubly closed-shell 
nuclei ^He and ^^O in the framework of the Goldstone expansion, starting 
from different realistic NN potentials. In order to rcnormalize their short- 
range repulsion, the high-momentum components of these potentials have 
been integrated out through the so-called Viow-k approach. We have em- 
ployed a criterion to map out the model space of the two-nucleon states in 
the HO basis according to the value of the cutoff momentum A [5] . 

To show the validity of this procedure, we have calculated the g.s. energy 
of ''He, with the CD-Bonn, N'^LO, and Bonn A potentials, comparing the 
results with the FY ones. We have found that the energy differences do 
not exceed 0.39 MeV. The limited size of the discrepancies evidences that 
our approach provides a reliable way to renormalize the NN potentials 
preserving not only the two-body but also the many-body physics. 

On the above grounds, we have performed similar calculations for ^^O 
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with the CD-Bonn and Bonn A NN potentials, and obtained converged 
results using model spaces not exceeding A^^max = 9. The rapid convergence 
of the results with the size of the HO model space makes it very interesting 
to study heavier systems employing our approach [17]. 
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